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Abstract 

We address the dynamics induced by collective atomic recoil in a Bose-Einstein condensate in 
presence of radiation losses and atomic decoherence. In particular, we focus on the linear regime 
of the lasing mechanism, and analyze the effects of losses and decoherence on the generation of 
entanglement. The dynamics is that of three bosons, two atomic modes interacting with a single- 
mode radiation field, coupled with a bath of oscillators. The resulting three-mode dissipative 
Master equation is solved analytically in terms of the Wigner function. We examine in details 
the two complementary limits of high-Q cavity and bad-cavity, the latter corresponding to the 
so-called superradiant regime, both in the quasi-classical and quantum regimes. We found that 
three-mode entanglement as well as two-mode atom-atom and atom-radiation entanglement is 
generally robust against losses and decoherence,thus making the present system a good candidate 
for the experimental observation of entanglement in condensate systems. In particular, steady-state 
entanglement may be obtained both between atoms with opposite momenta and between atoms 
and photons. 

PACS numbers: 42.50.Fx, 03.75.Gg, 42.50.Vk, 42.50.Dv, 03.67.Mn 
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I. INTRODUCTION 



The experimental realization of Bose-Einstein condensation opened the possibility to 
generate macroscopic atomic fields whose quantum statistical properties can in principle be 
manipulated and controlled [lj]. The system considered here to this purpose is an elongated 
Bose-Einstein Condensate (BEC) driven by a far off-resonant pump laser of wave vector 
k p = ujp/c along the condensate long axis and coupled to a single mode in an optical ring 
cavity. The mechanism at the basis of this kind of physics is the so-called Collective Atomic 
Recoil La 8 ie g (CARL)Q in Ki 8 full qUantiz ed version Q Q Q. In CARL the 8 e att e red 
radiation mode and the atomic momentum side modes become macroscopically occupied 
via a collective instability. A peculiar aspect of the quantum regime is the possibility of 
populating single momentum modes separated by Ap = 2hk p off the condensate ground 
state with zero initial momentum. The experimental observation of CARL in a BEC has 
been until now realized in the so-called superradiant regime ((J ?1 Q , i- e. without the optical 
cavity. In this case the radiation is emitted along the 'end-fire modes' of the condensate 
with very large radiation losses (in the mean field model, with n ~ c/L, where k is the 
cavity decay rate and L is the condensate length). In a recent work [5] it has been shown 
that atom-atom and atom-photon entanglement can be produced in the linear regime of 
CARL, in which the ground state of the condensate remains approximately undepleted. In 
this regime the atomic multi-mode system can be described by only two momentum side 
modes, with p = ±2Hk v . This source of entanglement has been also proposed for a quantum 

n n 

teleportation scheme among atoms and photons 10]. The results presented in |5| refer to 
the ideal case of a perfect optical cavity and an atomic system free of decoherence. However, 
in view of an experimental observation of entanglement, a detailed analysis of the sources 
of noise is in order, which in turn may be a serious limitation for entanglement in CARL 
[ll|. Also, it has not yet been proved that BEC superradiance experiments may generate 
entangled atom-photon states, as suggested in 0|. This issue is investigated for the first 
time in this paper, where we demonstrate the entangled properties of the atom-atom and 
atom-photon pairs produced in the linear stage of the superradiant CARL regime in a BEC. 

The aim of the present work is to analyze systematically, by solving the three-mode Master 
equation in the Wigner representation, the effects of losses and decoherence on the generation 
of entanglement. We will first investigate the effects of either a small atomic decoherence or 
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a finite mirror transmission of the optical cavity, and then analyze in details the generation 
of entanglement in the superradiant regime, where the cavity losses are important. 

The paper is structured as follows. In Section |H] we briefly review the ideal dynamics and 
derive the general solution of the Master equation. In Section UTTl we consider the evolution 
of the system starting from the vacuum and calculate the relevant expectation values, such 
as average and variance of the occupation number and two-mode squeezing parameters. In 
Section HVl the different working regimes are introduced and the dynamics analyzed, whereas 
in Section |3 we investigate three- and two-mode entanglement properties of the system as 
a function of loss and decoherence parameters. Section ED closes the paper with some 
concluding remarks. 



II. DISSIPATIVE MASTER EQUATION 

We consider a ID geometry in which a off-resonant laser pulse, with Rabi frequency 
Qq = dEo/h (where d is the dipole matrix element and Eq is the electric field amplitude) 
and detuned from the atomic resonance by A — UJp — Cl>q, is injected in a ring cavity 
aligned with the symmetry z-axis of an elongated BEC. The dimensionless position and 
momentum of the atom along the axis z are 9 = 2k p z and p = p z /2hk p . The interaction 
time is r = poo r t, where uj t = 2%k 2 /m is the recoil frequency, m is the atomic mass, p = 
(fio/2Ao) 2 ^ 3 (uj p d 2 N/Vheouj 2 ) 1 ^ 3 is the CARL parameter, iV is the number of atoms in the 
cavity mode volume V and €q is the permittivity of the free space. 

In a second quantized model for CARL jj, 3] the atomic field operator &(9) obeys the 
bosonic equal-time commutation relations &(9')] = 5(9-9'), = and the 

normalization condition is J 27r d9^f (9)^(9) = N. We assume that the atoms are delocalized 
inside the condensate and that, at zero temperature, the momentum uncertainty a Vz ~ h/a z 
can be neglected with respect to 2hk p . This approximation is valid for L ^> X p , where 
L is the condensate length and X p = 27c/k p is the laser radiation wavelength. In this 
limit, we can introduce creation and annihilation operators for an atom with a definite 
momentum p, i.e. ^{9) = J2 m c m (9\m) , where p\m) = m\m) (with m = — oo,...,oo), 
(9\m) = (l/\/27r) exp(im#) and c m are bosonic operators obeying the commutation relations 



3 



151 



[c m , c}J = and [c m , c n ] =0. The Hamiltonian in this case is 

H= £ {— & + iJ^(a ] cl l c n+1 - h.c.) j - 5a f a (1) 

n=— oo IP J 

where a is the annihilation operator (with [a, a'] = 1) for the cavity mode (propagating 
along the positive direction of the z-axis) with frequency uo s and 5 = {u p — u s ) / puo r is the 
detuning with respect to the pump frequency ui p . Let us now consider the equilibrium state 
with no photons and all the atoms at rest, i.e. with |^ ) = VN\0}- Linearizing around 
this equilibrium state and defining the operators a\ = c_ie 4<5r , a 2 = C\e~ %&T and a 3 = ae~ tSr , 
the Hamiltonian (JJJ reduces to that for three parametrically coupled harmonic oscillator 
operators: 

H = 8 + a\a2 — 8-a\ai + (^1 + 02)0-3 — (Si + 0-2)0.3 , (2) 

where 5± — 5 ± 1/ p. In Ref . H we have explicitly evaluated the state evolved from the 
vacuum of the three modes, |0i, 2 , 3 ), as 

where [hi) = (a|oj) with i — 1, 2, 3 are the expectation values of the occupation numbers of 
the three modes, related by the constant of motion C = fit — fi2 — n^. In this paper we extend 
our previous analysis to include the effects of atomic decoherence and cavity radiation losses. 
In this case the dynamics of the system in described by the following Master equation: 

-i H,g + 27iL[di]f) + 27 2 L[d 2 ]f? + 2kL[cl 3 ]q, (4) 



dg 
dr 



where 71, 72 and k are the damping rates for the modes Oj and L[ai] is the Lindblad super- 
operator 

L[di]g = diga\ - -o]a^ - -ga\di. (5) 

The atomic decay stems from coherence loss between the undepleted ground state with 
p z = and the side modes with p z = ±2hk v . In general, we assume that the two atomic 
modes may have different decoherence rates, depending on the direction of recoil [8(. The 
radiation decay constant is k = cT/2£, where T is the transmission of the cavity and £ is the 
cavity length. Through a standard procedure [l2j| . the Master equation can be transformed 
into a Fokker-Planck equation for the Wigner function of the state g, 

3 d% 
i=i 



W(ai, a 2 , «3, t) = /n^r e^-<^ X (ei,6,e3,T) , (6) 



where otj and 6 are complex numbers and x is the characteristic function defined as 



x(6,6,6 



Tr 



(7) 



where Dj(£j) = exp(6aj — is a displacement operator for the j-th mode. Using the 
differential representation of the Lindblad superoperator, the Fokker-Planck equation is: 

dW 



dr 



where 



u 



T 



u' T Au + c.c.j W + u' T Bu'*W 

d d d 



u 



IT 



\da\ da 2 da 3 
and A and D are the following drift and diffusion matrices: 



(8) 
(9) 



A 







p/2 



' 71 + ib- 

72 + i5+ 



^71 N 



D 



/ 



72 
0k 



(10) 



The solution of the Fokker-Planck equation (jHj) reads as follows 

W{u, t) = J d 2 u W{u , 0) G{u, r; uo, 0) (11) 

where W(uo, 0) is the Wigner function for the initial state and the Green function 
G(u, t; uo, 0) is the solution of Eq.© for the initial condition G(u, 0; u , 0) = 5^(u — u ). 
The calculation of the Green function, solution of Eq.fjHJ), is reported in detail in Appendix 

A and yields the following result: 
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G(u,t; u ,0) 



7r 3 det Q(r] 



exp {- [u - M(r)u ] f Q-\r) [u - M(r)u ]} . (12) 



where 



M( 



and 



Q(r) 



' fii{r) / 12 (r) / 13 (r) N 

-/i2(r) / 2 2(r) /2s(r) 
V /i3(t) -/ 23 (t) /3s(r) y 

dr' M(r') D M^(r') . 



(13) 



(14) 



In Eq. (jl3|) the complex functions fy, given explicitly in Appendix B, are the sum of three 
terms proportional to e tUkT , where ouj., with k = 1,2,3, are the three roots of the cubic 
equation: 



[u-S-i(k- 7+)] 
and 7± = (71 ± 7 2 )/2. 



'1 

- + i'J- 



+ 1 + zp7_ 







(15) 



III. EVOLUTION FROM VACUUM AND EXPECTATION VALUES 



Let now assume that the initial state is the vacuum. The characteristic function and the 
Wigner function at r = are given by 

=exp{-tfCo£} W(u,0) = fy* exp {-vtC^u} . , (16) 

where £ = (^1^2,^3) an d the covariance matrix is multiple of the identity matrix Co = |l. 
Since the initial state is Gaussian and the convolution in (|TT|) maintains this character we 
have that the Wigner function is Gaussian at any time r. After some algebra, we found that 
the covariance matrix is given by 

C(r) = Q(r) + ^M(r)Mt(r), (17) 

where the explicit form of the elements Cy = ((ui — (u)i)(uj — (u)j)*) in terms of the 
functions is reported in appendix |Bj Since the state is Gaussian, from (JT7J) it is possible 
to derive all the expectation values for the three modes. In particular, Ca = 1/2 + {hi), 
C12 = (^1^2)) C13 = an d C23 = (cl2&\)- The number variances and the equal-time 

correlation functions for the mode numbers are calculated from the forth-order covariance 
matrix G-ijki = ((ui — (u)i)(uj — (u)j)(uk — (u)k)*(ui — (u)i)*), which in turn is related to 
covariance matrix as follows: 

Gyki = CkiCij + CuCkj ■ (18) 

In particular, we have 

G m = (n?) + (hi) + \ (19) 

111 

Gijij = (hihj) + -(hi) + ~(hj) + - (i/ 3) (20) 

From Eqs.(|IEJ-<|2nj it follows that: 

v 2 (rii) = (n,)((n,) + l) (21) 



(2) _ (a|«l«A) ,,-yy, 

9i - { ^ )2 -I (-) 



(2) Y^i^j ) I C'ij I / \ 

9i ' j = WW) = 1 + WW) ' ( "' J)) 



where a 2 (rii) = (nf) — (hi) 2 , with i = 1,2,3, and i 7^ j in Eq.(J2HJ). The two-mode number 
squeezing parameter is calculated as [13]: 

= a 2 (hj -hj) = a 2 (n t ) +a 2 (n j ) -2\C l3 \ 2 

w (h i ) + (h j ) (h t ) + (h 3 ) ■ 1 ] 



We observe, from Eqs. (j21D and ()22|) that the statistics is that of a chaotic (i.e. thermal) 
state, as obtained in Ref. [5J] for the lossless case. If the two modes are perfectly number- 
squeezed, then £jj = 0, whereas if they are independent and coherent, ^ = 1. As it will be 
clear in the following sections, it is also worth to introduce also the atomic density operator 
for the linearized matter-wave field « [\/N + aie~ i{ - e+5T "> + a 2 e i( - e+ST ' ) ]/ V2rr, defined as 

n{9) = ¥(6)^(6) « — (l + Be- i{e+5r) + 5V (9+<5r) ) , (25) 

where B = (a\ + a 2 ) / \/N is the bunching operator, with (B) = and 

(&B) = hc u + C 22 + C l2 + C 21 ). (26) 



IV. ANALYSIS OF WORKING REGIMES 



We now investigate the different regimes of operation of CARL. For sake of simplicity, 
we will discuss only the case with 7i = 72 = 7, so that 7+ = 7 and 7_ = 0. In this case the 
cubic equation (fTB]) becomes: 

[ u - 5 - i (« _ 7 )] ^ - + 1 = (27) 

We will discuss two pairs different regimes of CARL, as defined in ref.^j], i.e.: i) semi- 
classical good-cavity regime (p> 1 and k <C 1); ii) quantum good-cavity regime (k 2 <C p < 
1); iii) semi-classical superradiant regime (p 3> \/2k > 1); iv) quantum superradiant regime 
(k 2 ^> \/2k > p). Also, we note that the case 7 = k worth a special attention. In fact, in 
this case Eq. (j27|) is independent on losses: the effect of decoherence is only a overall factor 
exp(— 77-) multiplying the functions fy, elements of the matrix M. Hence, it is expected 
that the case 7 = k will have statistical properties similar to those of the ideal case without 
losses, as it will be discussed below. 



A. CARL instability 

First, we investigate the effect of decoherence and cavity losses on the CARL instability 
in the different regimes. For large values of r the functions f\j of Eq. (jl3|) grow as exp(<?T), 
where g = — Imu — 7 is the exponential gain and u is the unstable root of Eq. (}2Tj) . with 
negative imaginary part. In figQwe plot g vs. 5 in the semi-classical regime (e.g. p = 100) 
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for the good-cavity case (k = 0) and 7 = 0.5,1,2 (figH^,), whereas the transition to the 
superradiant regime is shown in figCJ) for k — 1, 5, 10 and 7 = 0. The dashed line in figQ 
shows the gain for the ideal case k = 7 = 0. 

A similar behavior is obtained in the quantum regime shown in fig|21 where g is plotted 
vs. 5 for p = 0.2, k = and 7 = 0.2, 0.5, 1 (fig|2K quantum good-cavity regime) and for 
p = 1, 7 = and n = 0.5, 1,5 (figJ2b, quantum superradiant regime). Note that, unlike in 
the semi-classical regime, in the quantum regime the gain is symmetric around the resonance 
S = l/p (i.e. uj s = uj p — uj r ). Notice that in the case 7 = k, g = g^ — 7, where g^ is 
shown by a dashed lines in fig^andEl Whereas increasing k or 7 g tends to zero remaining 
positive for some value of S instead in the case 7 = k we have a threshold for g(°> = 7. 

B. Average populations and number squeezing parameter 

Figures El and 0] show the effect of losses, in the semi-classical regime, on the atomic 
population (hi), (a), and on the number squeezing parameter £i j2 , (b), plotted as a function 
of S for p = 100 and r = 2. Fig|3]shows the effect of the atomic decoherence on the the high- 
Q cavity regime (k = 0) for 7 = (dashed line), 0.5 and 1. We observe that increasing 7 the 
population of the mode 1 decreases and the number squeezing parameter £ 12 increases in the 
region of detuning where is less than one, i.e. where atom-atom number squeezing occurs. A 
similar behavior can be observed increasing the radiation losses in the semi-classical regime, 
as shown in figlH where (hi), (a), and £i )2 , (b), are plotted vs. 5 for 7 = 0, k = 0, 1, 5 
and p = 100. In both the cases, in order to observe number squeezing in the semi-classical 
regime, it is necessary to detune the probe field from resonance, as it was already pointed 
out in ref. p|. The inclusion of losses allows also to reach a steady-state regime when the gain 
g is negative. In this case, the covariance matrix C(oo) = Q(oo) becomes asymptotically 
constant. An example of this behavior is shown in figEl where (hi), (a), and £i i2 , (b), are 
plotted vs. r for p = 100 and 5 = 3.5. The dashed line shows the ideal case 7 = k = 0: 
because g^ = (as it can be observed from figHJ), the solution is oscillating and the two 
atomic modes 1 and 2 are periodically number squeezed. The dotted line of figH3 shows the 
case with k = and 7 = 0.2. Here, g = 0.025 and both the average population and the 
number squeezing parameter grow in time. Finally, the continuous line of figEl shows the 
case 7 = k = 0.5: the gain is g = —0.5 and the system reaches a stationary state in which 
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£1,2 = 0-7- This case is of some interest because a steady-state atom-atom number squeezed 
state is obtained in a linear system. 

Let now consider the effect of losses on the quantum regime. FigEl shows the average 
population {hi), (a), and the atom-photon number squeezing parameter £ 13 , (b), as a func- 
tion of r for 1/p = 5 = 5. Dashed lines in figEJa and b are for k = 7 = 0, the dotted lines 
are for k = and 7 = 0.15 and the continuous lines are for k = 7 = 0.15. We note that 
the atomic decoherence (i.e. 7) causes a drastic reduction of the number-squeezing between 
atoms and photons. However, choosing 7 = k < (where < yJp/2), we may keep 
£l ( 3 constant and less than one for a relatively long time, like in the ideal case case without 
losses. Notice that in the quantum regime the below-threshold regime (i.e. g < 7) is not of 
interest because the average number of quanta generated in each modes remains less than 
one. 



C. Superradiant regime 

In this section we present analytical results for the superradiant regime in the asymptotic 
limit |Ima;i|r ^> 1, where uj\ is the unstable root of Eq. (|27|) with negative imaginary part. 
For k^> \ui\ and assuming for simplicity 7 = 0, one root of Eq.(|27j) can be discharged as it 
decays to zero as exp (—act) and the other two roots may be obtained solving the following 
quadratic equation: 

9 uj + S + in 1 . , 

^ + ( f + «)»-i^ -^ = - (28) 

From Eq.()28J) it is possible to calculate explicitly the unstable root and evaluate asymptot- 
ically the expressions of the function appearing in Eq. (jl3|) . From them, it is possible to 
evaluate the expectation values of the occupation numbers in the semi-classical and quantum 
regimes. 



1. Semi- classical limit of the superradiant regime 



For k 3/2 > 1 > v^/P an d 5 = 0, the solutions of Eq. (|2*Hj) are =« (1 =F i)/V2k and 
the average occupation numbers are: 

.2 



(ni) 



P 
16k 



1 + 



/2k 
P 



,<2/k)V*t 



(29) 
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(ns> 



16k 



(30) 
(31) 



We observe that (ni) ~ (722) and (7x3) ~ {2/ pn){n\) <C (tii), so that the number of emit- 
ted photons is much smaller than the number of atoms in the two motional states. The 
asymptotic expression of the expectation value (|26|) of the bunching parameter is 



v 1 AN 



'2k 



P 



3 (2/ K )V2 7 



(32) 



Assuming that (B^B) approaches a maximum value of the order of one, then the maximum 
average number of emitted photons is about pN/2n 2 , whereas the maximum fraction of 
atoms gaining a momentum 2hk p is about p 2 /An. 



2. Quantum limit of the superradiant regime 

For 1 > pj ^/k and 5 = 1/p, the solutions of Eq. 

the average occupation numbers are: 



<*2> 
<*S> 



p 



'2k, 



d {p/k)t 



are u 12 ~ 1/p =F *p/ (2k) and 



(33) 



d {p/k)t 



\2yfc) 
2k 2 



(34) 
(35) 



In this case, (712,3) <^ an d (712) ~ (p/2) 3 { n 3)'- the average number of emitted photons 
is much less than the average number of atoms scattering a photon from the pump to the 
probe. Furthermore, the number of atoms making the reverse process, i.e. scattering a 
photon from the probe to the pump, can be larger than the number of photons scattered 
into the probe mode if p > 2, as it occurs in the current experiment on BEC superradiance 
0, 0] . In this regime the asymptotic expression of the expectation value of the bunching 
parameter is: 



1 

N 



H 


V2^) 







3 (p/«0 7 



(36) 



so that (7x3) ~ (pN/2K 2 )(B Jf B) , as in the semi-classical limit. The only difference is that 
in the quantum regime the maximum of (B^B) is 1/2, so that the maximum number of 
scattered photon in the quantum limit is half of that obtained in the semi-classical limit. 
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V. ENTANGLEMENT AND SEPARABILITY 



In this Section we analyze the kind of entanglement that can be generated from our 
system. First, we establish notation and illustrate the separability criteria. We also apply 
the criteria to the state obtained in the ideal dynamics. Then, we address the effects of 
losses. We study both the separability properties of the tripartite state resulting from the 
evolution from the vacuum, as well as of the three two-mode states that are obtained by 
partial tracing over one of the modes. The basis of our analysis is that both the tripartite 
state and the partial traces are Gaussian states at any time. Therefore, we are able to 
fully characterize three-mode and two-mode entanglement as a function of the interaction 
parameters 3 Q • 

A. Three-mode entanglement 

Concerning entanglement properties, three-mode states may be classified as follows }l(| : 
Class 1 : fully inseparable states, i.e. not separable for any grouping of the modes; 

Class 2 : one-mode biseparable states, which are separable if two of the modes are grouped 
together, but inseparable with respect to the other groupings; 

Class 3 : two-mode biseparable states, which are separable with respect to two of the three 
possible bipartite groupings but inseparable with respect to the third; 

Class 4 '■ three-mode biseparable states, which are separable with respect to all three bipartite 
groupings, but cannot be written product state; 

Class 5 : fully separable states, which can be written as a three-mode product state. 

Separability properties are determined by the characteristic function. In order to simplify 
the analysis we rewrite the characteristic function (JJJ) in terms of the real variables x T = 
(xi, x 2 , x 3 , yi, y 2 , y 3 ) with ^ = 2~ 1 / 2 (?/ j - ixj), j = 1, 2, 3. We have 

X (x)=exp|-lx T Vx} , (37) 



where 

/ 

V = 2A 



A -B 
B A 



A , (38) 
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with A = Diag(— 1, 1, 1, 1, 1, 1) and 



A = 3?eC B = 3mC 



(39) 



and where we omitted the explicit time dependence of the matrices. The entanglement 
properties of the three-mode state are determined by the positivity of the matrices 



J = 1,2,3 




where A 1 = Diag(l, 1, 1, -1, 1, 1), A 2 = Diag(l, 1, 1, 1, -1, 1), A 3 = Diag(l, 1, 1, 1, 1, -1) and 
J is the symplectic block matrix 



(40) 



I being the 3x3 identity matrix. The positivity of the matrix Tj indicates that the j-th 
mode may be factorized from the other two. Therefore, we have that i) if Tj < Vj the 
state is in class 1; ii) if only one of the Tj is positive the state is in class 2; iii) is only two 
of the Tj are positive the state is in class 3; iv) if Tj > 0, Vj then the state is either in class 
4 or in class 5. 

The covariance matrix V can be written as 

In A V n ^ c \ 



V 



where 



Q 


-A —B V 


£ 


-A 


n c v o 


-T 


-B 


C 1ST 








V £ g A 


B 


V 


T AH 


C 


s 


-T B C 


1 


V 


= 2 Sra C12 


Q = 


£ 


= 2 3m C13 


H = 


T 


= 2 3m C23 


1 = 



(41) 



(42) 



and the matrix elements CV, are reported in Appendix [BJ Let us first consider the ideal case, 
when no losses are present. In this case we can prove analytically that the evolved state (JH1) 
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is fully inseparable. In fact, we have that 



A 



2 yj (n 2 )(l + (hi)) COS 02 



V 



2yJ(n 2 ){l + (ni))sin02 



2 ^ (n 3 )(l + (ni)) cos 3 



2 v /(n 3 )(l + (ni))sin0 3 



c 



2yJ(h 2 )(h 3 ) cos(0 3 - 02 ) 



2yJ(h 2 ){h 3 ) sin(0 3 - 2 ) , 



(43) 



from which, in turn, it is straightforward to prove that the minimum eigenvalues of the 
matrices Tj are always negative. In the non ideal case, when 7 or k are different from zero, 
the expressions given in Eqs. (|42jl and accordingly the minimum eigenvalues of matrices Tj 
should be calculated numerically. In Fig. [71 |H] and 01 the minimum eigenvalues of matrices 
r 1; T 2 and T 3 are plotted in the semi-classical regime, with p = 100. In this regime we 
can observe that modes 1 and 2 remain non separable from the three mode state even for 
large values of atomic decoherence 7 and radiation losses k. Instead inseparability of mode 
3 is not so robust especially in presence of some atomic decoherence. In Fig. El and ^2 
the minimum eigenvalues of matrices Y\ and T2 are plotted in the quantum regime, with 
p = 0.2. The minimum eigenvalue of the matrix T 3 is not reported in the figure since the 
behavior is similar to that of V\. In this regime we can observe that modes 1 and 3 remain 
non separable from the three mode state even for large values of atomic decoherence 7 and 
radiation losses k. Instead inseparability of mode 2 is very sensible especially in presence of 
some radiation losses. In any case in the quantum regime the three eigenvalues increasing 7 
and k approaches to zero but remain negative. In the semi-classical regime the eigenvalue 
of T 3 that corresponds to photonic mode 3 becomes positive increasing 7. 

B. Two- mode entanglement 

In experimental conditions where only two of the modes are available for investigations, 
the relevant piece of information is contained in the partial traces of the global three-mode 
state. Therefore, besides the study of three-mode entanglement it is also of interest to 
analyze the two-mode entanglement properties of partial traces. At first we notice that the 
Gaussian character of the state is preserved by the partial trace operation. Moreover, the 
covariance matrices V^- of three possible partial traces Qij = Tik [g] , i 7^ j 7^ k can be obtained 
from V by deleting the corresponding k-th and k + 3-th rows and columns. The Gaussian 
character of the partial traces also permits to check separability using the necessary and 
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sufficient conditions introduced in Ref. [15j, namely by the positivity of the matrices Sy 
and S^„- that are obtained by deleting the k-th and k + 3-th rows and columns either from 
Tj or Tj. Since they differ only for the sign of some off-diagonal elements it is easy to prove 
that they have the same eigenvalues. Therefore, we employ only Sy in checking separability. 
The matrices are given by 

( r. a a Ti \ 



'12 



'13 



'23 



V 



r> 
y 


A 


i 




A 

—A. 


1-1 
it 


— u 


i 


— % 


— L> 


n. 

y 


A 

— J\ 


V 


—i 


-A 


H 


Q 


-B 


i 


£ \ 


-B 


C 


-8 


i 


—i 


-8 


Q 


-B 


8 


—i 


-B 


1 J 


n 


C 


i 


-T 


c 


£ 


-C 


% 


—i 


-C 
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-T 
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(44) 



(45) 



(46) 



In ideal conditions with 7 = k = the minimum eigenvalues of Si^, k = 2, 3 are given by 

Vik = (ni) + (n k ) - ^4,{h k ) + ((ni) + (h k )) 2 (47) 

and thus are always negative. On the contrary, the minimum eigenvalue of S23 is given by 

V23 = 1 + (ni) + v /(l + (n 1 )) 2 -4(n fc ) (48) 

where {hk) = max((n2), (^3)). Note that 7723 is always positive. Therefore, after partial 
tracing we may have atom-atom entanglement (entanglement between mode ai and mode 
02) or scattered atom-radiation entanglement (entanglement between mode a± and mode 03) 
but no entanglement between mode and mode a 3 . 

For r> 1 we know the asymptotic expressions for populations in the ideal case without 
losses j^J, so we can obtain the stationary value of 7]^ as 

2{h k ) 



{hi) + {h k ) 
14 



(49) 



In the high-gain semi-classical regime (p 3> 



1) S, 



(hi) 
(n 2 > 
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so that 
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In the high-gain quantum regime (p < 1), 
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so that 
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Vl3 
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16 + p 3 



(50) 

(51) 
(52) 

(53) 

(54) 

(55) 
(56) 

(57) 



4 + p 3 

In the non ideal case, when 7 or k are different from zero, the minimum eigenvalues of 
matrices Si 2 and S13 can be easily obtained numerically. In Fig. ^] and El the minimum 
eigenvalues of matrices S12, S13 are plotted for the semi-classical regime. We can observe that 
the atom-atom entanglement of the reduced state 12 is robust, as the minimum eigenvalue 
remain negative increasing atomic decoherence 7 and radiation losses k. On the contrary 
atom-photon entanglement of the reduced state 13 is more sensitive to noise: the eigenvalue 
remains negative increasing k and become positive in presence of some atomic decoherence. 

In Fig. Eland^Jare plotted the minimum eigenvalues of matrices S12, S13 in the quantum 
regime, with p = 0.2. Here the atom-photon entanglement in the state 13 is robust while 
atom-atom entanglement of the state 12 is not. The minimum eigenvalue always remains 
negative, but it starts from a very small absolute value and approaches very fast to zero 
increasing k and 7. 



VI. CONCLUSIONS 



We have investigated how cavity radiation losses and atomic decoherence influence the 
generation of two (atom-atom or atom-radiation) and three mode entanglement in the collec- 
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tive atomic recoil lasing (CARL) by a Bose-Einstein condensate driven by a far off-resonant 
pump laser. The atoms back-scatter photons from the pump to a weak radiation mode 
circulating in a ring cavity, recoiling with opposite momentum ±2hk p along the ring cav- 
ity axis. Our analysis has been focused to the linear regime, in which the ground state of 
the condensate remains approximately undepleted and the dynamics is described by three 
parametrically coupled boson operators, corresponding to the radiation mode and two con- 
densates with momentum displace d by ± 2hk p . The problem resembles that of three optical 
modes generated in a medium |17[ and thus our results may have a more general interest 
also behind the physics of the BEC. We have solved analytically the dissipative Master equa- 
tion in terms of the Wigner function and we have investigated the entanglement properties 
of the evolved state. We found that three-mode entanglement as well two-mode atom-atom 
and atom-photon entanglement is generally robust against cavity losses and decoherence. 
The analysis has been focused of the different dynamical regimes, the high-Q cavity regime, 
with low cavity losses, and the superradiant regime in the so-called 'bad-cavity limit'. We 
have found that entanglement in the high-Q cavity regime is generally robust against ei- 
ther cavity or decoherence losses. On the contrary, losses seriously limit atom-atom and 
atom-radiation number squeezing produ c tlon in CARL Q. Conce^ the ^ 
regime, atom-atom entanglement in the semi-classical limit is generally more robust than 
atom-radiation entanglement in the quantum- limit. Finally, we have proved that the state 
generated in the ideal case without losses is fully inseparable. We conclude that the present 
system is a good candidate for the experimental observation of entanglement in condensate 
systems since, in particular, steady-state entanglement may be obtained both between atoms 
with opposite momenta and between atoms and photons. 
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APPENDIX A: SOLUTION OF THE FOKKER-PLANCK EQUATION 



In order to solve Eq.fjHJ) for the Green function G(u, t; Uq, 0) it is helpful to first perform 
a similarity transformation to diagonalize the drift matrix A: 



A = SAS 



-i 



diag{AiA 2 A 3 }, 



(Al) 



where the complex eigenvalues Xj of A (with j = 1,2,3) are obtained from the characteristic 
equation 

det(A-AI) = 0, (A2) 

I is the 3x3 identity matrix and the columns of S _1 are the right eigenvectors of A with 
det(S) = 1. Solving Eq. (jA2|) we obtain Xj = i(uj — 5) — 7+, where Uj are the three roots of 
the cubic equation (|T5j) . whereas the eigenvectors of A corresponding to the j-th eigenvalue 
are 



where (3 = 1/p + «7_ and 



l 



M 2 



l 



M% 



UJ 2 — 0J 3 W\— UJ 3 OJl — UJ 2 

Explicitly calculating the inverse matrix of S~ l we have 

/ (a 22 a 23 /M 1 ) (a 12 a 13 /M) -M 2 M 3 /M^ 



(A3) 



(A4) 



-i^fp/2 (a 23 a 21 /M 2 ) -iJp/2{a u a 13 /M 2 ) M 1 M 3 /M 2 



(A5) 



where a 



Kl 



\ iyfp/2 (a 21 a 22 /M 3 ) (a 12ail /M 3 ) -M X M 2 /Mz, ) 

(dLj)i. Now we transform the Fokker-Plank equation (JSJ) in the new variable 



v = Su. From (jAl|) we obtain 

u /T Au = u ,T (S- 1 AS)u = v' T Av 

jTa\nfaT„j\* iT-f\.j* 



(A6) 

v ,T S)D(S r v')* = v ,T BV*, (A7) 

where D = SDSt, S f = (S r )* and v ,T = u^S" 1 . Using (jMI) and Eq. © becomes 

dW 



Or 



- lv' T Av + c.c. ) W + V 1 By'*W, 



(A8) 



where W(v,t) = W(S 1 v,r). Eq. ()A8j) is a linear Fokker-Plank equation with diagonal 
drift. Introducing the Fourier transform 



' / " k IU(v)exp(k* T v-k T v*) ) 



7T 



(A9) 
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Eq. (|A8|) becomes 



where 



^ = (k* T Ak'* + k T A*k') U - (k* T Dk) U. 

d d d \ 



(fa, k 2 , fa) 



IT 



(A10) 



(All) 



V dki dk 2 dfa J 

The Fourier transform of initial condition of the Green function G(v, 0; Su , 0) = 5 3 (v — Su ) 
is 

U(k, 0) = exp [k* T Su - k T (Su )*] . (A12) 

Eq. flA10|) is now solved using the method of the characteristics. Since A is diagonal the 
subsidiary equations are 



dr dk\ 



dkl 



dk% 



dU 



and have solutions 



1 -X 1 k* 2 -\ 2 k* 2 -X 3 k* 3 (-k* T Dk)[7 



k = e A r c = constant. 



Then 
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:* T (e" Ar De 
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(A13) 



(A14) 



(A15) 



where (By) denotes the matrix with elements By, and we find, using Eq. (|A14|) . 



In U = k 



D 



Ai + A* 



I _ e (Ai+A;)r 



k + constant 



It follows that 



U exp |k* T Qk| = constant 
where Q is the 3x3 matrix with elements 

A, 



\i + A* 



1 - e 



(Ai+A*)r 



dT'biJ^W. 



(A16) 



(A17) 



(A18) 



Thus, from Eqs. ()A14|) and (|A17)1 . the solution for U takes the general form 



U(k, r) = $(e A r k) exp{-k* J Qk} 



(A19) 



where $ is an arbitrary function. Choosing $ to match the initial condition (jA12|) . we find 
U(k, t) = exp {k* T (Se Ar u ) - k T (Se Ar u )*} exp {-k* T Qk} . (A20) 
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In the argument of the first exponential on the right-hand side we have used 
exp(Ar)S = Sexp(Ar). Inverting the Fourier transform we obtain 

G(v, r; Su , 0) = 1 - exp { (v - Se A ^u ) j Q 1 (v - Se Ar u ) j 
tt 6 det Q L ' > 

and so transforming back the variables 

G{u, r; u , 0) = ^ 3 — - exp | (u - e Ar u ) Q 1 (u - e Ar u ) | . 



II) to write 



(A21) 



(A22) 



where 



Q = S _1 Q (S _lVf 



dr e D (e 



(A23) 



APPENDIX B: ELEMENTS OF THE MATRICES M, EQ. flUD, AND C, EQ. (fT7|> 

The expressions of the functions which appear as elements of the matrix M, Eq. (|13jh 



arc 
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3 ft 
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(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 



where a = 6 + i(n — 7+), /3 = 1/p + ry_, Aj = (u^- — uj k ){ujj — cj m ) (with j ^ k ^ m) and 
c^i, co> 2 and u 3 are the roots of the cubic Eq.([T5|). It is possible to show that /y(0) = <5y in 
order to satisfy the initial condition M(0) = I. 

The explicit components of the covariance matrix 



C(r) = Q(r) + -M(r)Mt(r), 



(B7) 
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where M and Q are denned in (|T3|) and (JHJ), are 
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with C tj = C£. 









, where f(lt) 

without losses. As shown in Ref. [5(, they satisfy the following relations: 



In the special case 7+ = 7 = k and 7_ = 0, = e yT fij , where /^ is the solution 
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I ill I ~ 1 

i/i ( 2 0) r+i 

/g> (/£>)• 

-A ( 2 0) (/l ( 3 0) )* 



(o), 2 + | /3 (o)| 2 



I/23 
Iil2 I 



+ i/i ( 3 0) r 



i/ 2 ( 2 0) r+i/2 ( 3 0) r 

-/l ( 2 0) (/2 ( 3 0) )* + /l ( 3 0) (/3 ( 3^ 

/^(/gV + ^aSV 

f(0)/ f (0K* , f (0)/ f (0) v 
J22 W23 I t i23 U33 J 



(B14) 
(B15) 
(B16) 
(B17) 
(B18) 
(B19) 



Using Eas. (lBT4l) - (lBT6jl in Ea.(lB8l- (lBT0jl and C u = 1/2+ (hi), we obtain that 

(fit) = (n 2 > + (n 3 ) 

and 



d(rij) d(nf^) 



-2jT 



(B20) 



(B21) 



dr dr 

where (nj '') are the expectation values of the occupation numbers of the three modes in the 
ideal case without losses. 
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FIG. 1: Growth rate g = — Imu; — 7 vs. 5 for the unstable root of the cubic equation (J15|) in the 
semi-classical limit, p = 100. In (a) k = and 7 = 0.5, 1,2; in (b) 7 = and k = 1, 5, 10. The 
dashed lines represent the case k = 7 = 0. 
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FIG. 2: Growth rate g = — Imu; — 7 vs. 5 for the unstable root of the cubic equation (J15|) in the 
quantum limit. In (a), p = 0.2, k = and 7 = 0.2,0.5,1; in (b), p = 1, 7 = and k = 0.5,1,5. 
The dashed lines represent the case n = 7 = 0. 
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FIG. 3: Semi-classical regime with k = 0: (ni), (a), and £1,2, (b), vs. 5 for p = 100, r = 2, 7 = 
(dashed line), 0.5 and 1. 
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FIG. 4: Semi-classical regime with 7 = 0: (ni), (a), and £1,2, (b), vs. 5 for p = 100, r = 2, re = 
(dashed line), 1 and 5. 
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FIG. 6: Quantum regime for l/p = 5 = 5: (ni), (a), and £1,3, (b), vs. r for 7 = 0, re = (dashed 
line), for 7 = 0.15, re = (dotted line) and for 7 = re = 0.15 (continuous line). 
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FIG. 7: Semi-classical regime for p = 100 and 5 = 0.01: minimum eigenvalue of matrix T\ for 
k = and different values of 7 (a) and for 7 = and different values of k (b). 
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FIG. 8: Semi-classical regime for p = 100 and 5 = 0.01: minimum eigenvalue of matrix T2 for 
k = and different values of 7 (a) and for 7 = and different values of k (b). 
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FIG. 9: Semi-classical regime for p = 100 and 5 = 0.01: minimum eigenvalue of matrix 1^ for 
k = and different values of 7 (a) and for 7 = and different values of k (b). 
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FIG. 10: Quantum regime for p = 0.2 and 6 = 5: minimum eigenvalue of matrix Ti (or T3, see the 
text) for k = and different values of 7 (a) and for 7 = and different values of k (b). 
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FIG. 11: Quantum regime for p = 0.2 and 5 = 5: minimum eigenvalue of matrix T2 for k = and 
different values of 7 (a) and for 7 = and different values of k (b). 






0.0j 




-0.2 




-0.4 


CM 




£T 


-0.6 




-0.8 




-1.0 









0.0 - 




-0.2 




-0.4 


CM 




-0.6 




-0.8 




-1 n 





K=0 


K=10 




K=100 







4 6 
X 



8 10 







II II II 


t 


(b) 







4, 6 



8 10 



FIG. 12: Semi-classical regime for p = 100 and 5 = 0: minimum eigenvalue of matrix Su for k = 
and different values of 7 (a) and for 7 = and different values of k (b). 
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FIG. 13: Semi-classical regime for p = 100 and 5 = 0: minimum eigenvalue of matrix Sis for k = 
and different values of 7 (a) and for 7 = and different values of k (b). 
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FIG. 14: Quantum regime for p = 0.2 and 5 = 5: minimum eigenvalue of matrix S12 for k = and 
different values of 7 (a) and for 7 = and different values of k (b). 
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FIG. 15: Quantum regime for p = 0.2 and 5 = 5: minimum eigenvalue of matrix S±s for k = and 
different values of 7 (a) and for 7 = and different values of k (b). 



30 



